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Abstract 

We study the statistics of local field distribution solved by the Green's- 

function formalism (GFF) [Y. Gu et al., Phys. Rev. B 59 12847 (1999)] in 

the disordered binary resonant composites. For a percolating network, the 

inverse participation ratios (IPR) with q = 2 are illustrated, as well as the 

typical local field distributions of localized and extended states. Numerical 

calculations indicate that for a definite fraction p the distribution function of 

IPR Pq has a scale invariant form. It is also shown the scaling behavior of 

the ensemble averaged (Pq) described by the fractal dimension Dq. To relate 

the eigenvectors correlations to resonance level statistics, the axial symmetry 

between D2 and the spectral compressibility x is obtained. 
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The statistics of eigenfunctions of the random matrix with orthogonal, unitary and sym- 
plectic symmetry has been well reviewed by Guhr^ and Mirlin^. A strong fluctuation of 
eigenfunctions at the critical point is one of the prominent hallmarks of Anderson metal- 
insulator transition. These fluctuations can be represented by a set of inverse participation 
ratios (IPR) Pq = J (i'^r|\E'(r)p^. Wegner put forward from the renormalization group treat- 
ment of the supersymmetry S model in 2 -|- e dimensions that IPR show the multifractal 
scahng^'^ with respect to the system size L, Pq — AL~^''^'^~^\ in a good metal, while in an 
insulator, Pg oc L°. For the critical power-law random bounded matrix (PRBM) ensemble^ 
with (3 = 1, not only the distribution function of the critical Pq was found to have a scale 
invariant form, but also the scaling exponent Dq was obtained from the shift of the distribu- 
tion of IPR with the size A^". At the criticahty, the multifractality of the ensemble- averaged 
(Pq) has been derived analytically in the limits of weak and strong couplings, as well as 
calculated numerically in the full range of couphngs^. Recently, local field distribution of 
resonant composites has attracted great interest. Various optical eigenmode localizations, 
such as surface-plasmon modes in the colloid clusters^, localized dipolar excitations on the 
roughly nanostructured surfaces^, and selective photomodification in the fractal aggregates 
of colloidal particles^, were reported to be relevant to the specific filed distribution. It was 
found that the local field at the percolation threshold fiuctuates from one optical eigenmode 
to another^°'^^. Some attempts to investigate the scaling of the successive local electric field 
moments have been made^^~^^. However, it is difficult to take into account the structure 
sensitivity to local field distribution when resonance happens. In order to study the statis- 
tics of resonant composites, one must answer the following questions. How do the optical 
eigenmodes distribute in the resonant area? For each mode, how about its local field distri- 
bution? And how to understand the fiuctuations and scaling of the local field distribution? 
It is known that the resonance spectrum and local field distribution for each eigenmode can 
be solved analytically using the Green 's-f unction formahsm (GFF)^^'^^. The first question 
has been answered in view of level spacing statistics^^'^^. In this work, local field distribu- 
tions of typical localized and extended eigenstates are illustrated, and to answer the third 
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question, the statistics of local field distribution is numerically studied from the random 
matrix theory (RMT). 

In this work, a binary network is considered where the impurity bonds with admittance 
ei are employed to replace the bonds in an otherwise homogeneous network of identical 
admittance 62- The admittance of each bond is generally complex and frequency-dependent. 
All the impurity bonds construct the clusters subspace. Using the GFF, the eigenvalues 
and eigenvectors of Green 's-matrix M can be solved analytically Because M maps 
the geometric configuration of the clusters subspace, it is a random matrix where RMT 
can set in. Its solutions summarize all the geometric resonances of the clusters subject to 
the external sources and in the quasistatic hmit. The element of M is defined as Mx,y = 
SzGC(y) (G^x,y " G^x,z) vo. which z G C (y) means that the jointing points z and y belong to 
the clusters subspace and are the nearest neighbors, and Gx,y is the Green's function of 
Laplace operator on the infinite square, i.e., — AG'x,y = 5x,y with Gx.x = 0. More clearly, 
Mx,y describes the interaction between x and y, and is closely related to the "environment" 
or the nearest neighbors of y. Different from the quantum systems, local field distribution is 
not directly equal to the eigenvector of M. When resonance happens, the local field can 
be expressed as a separate product of one sum of right eigenvectors of M and another sum of 
left eigenvectors that depends on the source term. Motivated by that local field distribution 
may have the general statistical features, this paper mainly devotes to the fiuctuations and 
scaling of the IPR. 

In the following, we start from solving Green's-matrix M numerically in the disordered 
binary networks. The distribution functions of the IPR are calculated for system size ranging 
from L = 16 to L = 38 and for various values of p. For the ensemble averaging over 1000 
samples are used in the percolating case at L = 30, for each of which there are more than 750 
nontrivial eigenstates to be produced. In the calculation of the scaling behavior of the IPR, 
{Pq) is averaged over 7 x 10^ right eigenvectors. The difference between using 1000 samples 
of 30 X 30 system and using a large system, say 300 x 300, is that there is a minor correction 
in the calculations. Numerical calculations indicate that the distribution function of IPR Pq 
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has a scale invariant form for a definite p within the interval [0.3,0.5]. It is also found the 
scaling behavior of the ensemble averaged (Pq) described by the fractal dimension Dq. The 
multifractal scaling of (Pq) with respect to L is obtained for p = 0.5 and p = 0.4, as well 
as the fractal behavior of (Pq) for p = 0.3. Note that here L represents the scale of square 
network, while in Refs. and |0, is the size of NxN matrix. To relate the level statistics 
to eigenvectors correlations, the symmetry between D2 and the spectral compressibility x is 
obtained. 

In order to have a direct observation, figure 1 displays the values of IPR of right eigenvec- 
tors with g = 2 for a 40 X 40 percolating network. The peaks represent the localized states 
and valleys the extended states. It is shown that peaks and valleys randomly distribute in 
the resonant area. The localized states incline to accumulate near s ^ and s 1.0, while 
at about s 0.5, more extended states are found. Various optical excitations^^^, optical 
nonlinear enhancements^^'^^, and eigenmode localizations^^ are typical localized states with 
the high values of IPR. For a localized state, it is found that the value of P2 of right eigen- 
vector is always high, as well as the IPR of potentials and electric fields^^. Therefore, it is 
reasonable to represent the localization or extension of eigenstates only by the IPR of right 
eigenvectors. When resonance happens, the inhomogeneous local field around the impurity 
metallic clusters leads to a large enhancement in the effective linear and nonlinear optical 
responses. It is expected that the high values of IPR, or the localized states, correspond 
to the strong nonlinear optical enhancements. There are an anomalous absorption in the 
infrared area in percolating metal-dielectric thin films^^'^° and a large optical nonlinear en- 
hancement in the region of high frequencies^^ . In ref . , the separation of the absorption 
peak from nonlinear enhancement peak is discussed. Different from the causes in the dilute 
anisotropic networks^^, it originates in the self-similar percolating structure. 

Then, corresponding to the different values of IPR in Fig. 1, local field distributions 
of typical localized and extended states are illustrated in the forms of 3D plots. Fig. 2(a) 
shows a very localized state with the high value of P2- It is seen that when one eigenmode 
is excited, the rest of the volume remains almost unexcited^. For simplicity, those strong 



localized fields are called "hot spots". In this figure, the brightest "hot spot" is found. In 
Fig. 2(b), an extended state with the low value of P2 is plotted, where local fields are not 
uniform and many smaller and weaker "hot spots" are found. In Ref . [14] , we also shown an 
intermediate state with the average value of P2, where we can see more "hot spots" and local 
fields are very inhomogeneous. "Hot spots" are found to be sensitive to the admittance ratio 
h{= ti/t2) and external fields. They also lead to the photomodification^ and huge enhanced 
nonlinear responses^^'^^. 

When we plot the IPR of several samples with different size, it is intuitive that the 
general feature of the distribution function of the IPR should exist. Figure 3 displays the 
distributions of ln(Pq/ {Pq)) for p = Pc = 0.5, where Pg is normalized by {Pg) and q is ranged 
from 2 to 8. It is seen that the shape of the distribution function changes with increasing q. 
It is also shown that for the distribution function of ln(Pg/ (Pg)) the change of p leads to the 
shift and deformation of the shape^'^. This is different from the recent argument of Mirlin 
in which different parameter b leads to the changes of the distribution function of In(Pg) for 
the critical PRBM ensemble^. 

Figure 4 displays the distributions of ln(P2) with the system size ranged from 16 to 32 
stepped by 4 ioi p = Pc = 0.5, i.e., the percolating case. While for p = 0.4 and p = 0.3, the 
distribution functions of ln(P2) with the different system size are also calculated^^. A scale 
invariant form of the IPR distribution is found for the critical case Pc = 0.5, as well as for 
the noncritical p = 0.4 and p = 0.3. Note that in the Refs. and 0, only at the criticality 
of PRBM ensemble, a scale invariant form of the distribution of the IPR is investigated. 
The shift of the curves implies a scaling of distribution function of the IPR with respect to 
L. There exist different scaling exponents Dg for the different values of g at p = 0.5, 0.4, 



and 0.3, as verified by the ensemble averaging (Pg) in ref. [Q. All curves are not strictly 
overlapped when we shift them by In (Pg), especially for the large q. One possible reason is 
that the sample is not large enough. Therefore, from Figures 3 and 4 we can conclude that 
for a definite p the fiuctuations of the IPR have the general statistical features. 

Different from the eigenf unctions of the Anderson metal-insulator Hamiltonian and 



PRBM ensemble, the right eigenvectors of Green's-matrix M are the main product of the 



local field V for the binary composites. In ref. we found the scaling of {Pq) with re- 
spect to the size L for pc = 0.5, p = 0.4 and p = 0.3. By the relation, {Pq) = L""^', the 
scaling exponent Tg can be calculated through the slope of In (Pg) with ln(L). The results 
are shown in Fig. 5. For p = 0.5 and p = 0.4, we observe the multifractal scaling of (Pq). 
Note that for p = 0.4 it is not a critical case. While in Ref. |P and Ref. the multifractal 
scaling is investigated only in the critical ensembles. However, for p = 0.3, only the fractal 
behavior of the ensemble averaged (Pq) is found and for q > 2 the extended line of Tq can be 
approximately approaching the original point (0, 0). By the relation Tq = Dq x (g — 1), the 
values of Dq are obtained. Here D2 is not closely related to the spatial geometry as claimed 
in Ref. |2ll . 



Finally, to connect with the level spacing statistics, the scaling exponent D2 and spectral 
compressibility^^ X are shown in Fig. 6, and the axial symmetry according to their mean 
values is found when p is ranged in the interval [0.3,0.5]. The result apparently disagrees 
with the claim of Ref. where D2 and x satisfy D2 + 2X = 1 for the large b, and 2D2+X = 1 
for the small b. It is also violated with the recent statement^^, where between D2 and x the 
exact relation holds, x = {d — D)/2d. 

In a summary, we have not only presented a detailed study on the local field distribution 
near resonance, but also the statistics of local field distribution by means of the IPR for the 
disordered binary networks. On the qualitative level, the main findings are summarized as 
follows: 1) For the percolating composites, the IPR of the right eigenvectors are illustrated, 
as well as the typical local field distributions of the localized and extended states in the 
forms of 3D plots. When resonance happens, the fields are localized within the impurity 
metallic clusters. For a percolating network, a large enhancement in the effective linear and 
nonlinear optical responses is also found^^. 2) The distribution functions of the IPR Pq have 
a scale invariant form in the limit of the large system size L for the critical case, pc = 0.5, 
as well as for the noncritical case, p = 0.4 or p = 0.3. 3) The multifractal scaling of the 
ensemble averaged (Pq) with the system size L is obtained for pc = 0.5 and p = 0.4, while 



for p — 0.3, only the fractal behavior is found. 4) To relate the eigenvectors correlations 
to level spacing statistics, the axial symmetry between D2 and x is obtained. Finally, it 
is worthwhile to emphasize on the extent of the universality of the IPR distribution at the 
noncritical case, such as p = 0.4 or p = 0.3. In the previous works^'^, only at the criticality, 
the fluctuations and scaling of the IPR were investigated. 
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FIGURES 

FIG. 1. IPR of right eigenvectors with q = 2 aX the percolating threshold Pc- Here we use a 

40 X 40 sample with 1331 resonances. 

FIG. 2. The typical local field distributions of localized and extended eigenstates. a) with a 
high value of IPR; b) with a low value of IPR. 

FIG. 3. Distribution functions of ln(Pq/ {Pq}) with respect to ln(Pq/ {Pq}) for q ranged from 
2 to 8 and for the percolating networks. Here L = 32 and the distributions are averaged over 900 
samples. 

FIG. 4. Scale invariant form of distribution functions ln(P2) with respect to ln(P2) for pc = 0.5. 

FIG. 5. Diagram of multifractal scaling Tq for p = 0.5, 0.4, and 0.3. For p = 0.3, the extended 
line of Tq can be approaching to the point (0,0) when q > 2, therefore only fractal behavior is 
found. 

FIG. 6. The fractal dimension D2 (circles) and the spectral compressibility x(squares) as a 
function of the parameter p. 
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